Memory Effects in Granular Material 
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We present a combined experimental and theoretical study of memory effects in vibration-induced 
compaction of granular materials. In particular, the response of the system to an abrupt change 
in shaking intensity is measured. At short times after the perturbation a granular analog of aging 
in glasses is observed. Using a simple two-state model, we are able to explain this short-time 
response. We also discuss the possibility for the system to obey an approximate pseudo-fiuctuation- 
dissipation theorem relationship and relate our work to earlier experimental and theoretical studies 
of the problem. 
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Granular materials comprise an important class of 
complex systems whose simple fundamental mechanics 
gives rise to rich macroscopic phenomenology (lj. Re- 
cent experiments on granular compaction [g,|3| suggest 
they are an ideal system for studying jamming, a phe- 
nomenon lying outside the domain of conventional statis- 
tical physics, yet highly reminiscent of glassiness. These 
studies showed that a loose packing of glass beads sub- 
jected to vertical "tapping" slowly compacts, asymptot- 
ing to a higher steady state packing fraction. This "equi- 
librium" packing fraction is somewhat lower than the 
random close packing limit, p rcp « 0.64, and is a decreas- 
ing function of the vibration intensity, typically parame- 
terized by r, the peak applied acceleration normalized by 
gravity, g. The relaxation dynamics are extremely slow, 
taking many thousands of taps for the packing fraction, 
p, to approach its steady state value. During this evolu- 
tion, p increases logarithmically with the number of taps, 
t, which is typical for self-inhibiting processes Q. The 
average time scale r of the relaxation decreases with T, 
and in this sense the shaking intensity plays, at least qual- 
itatively, the role of temperature. For small T, the relax- 
ation rate becomes so slow that the system cannot reach 
the steady state density within the experimental time 
scale. It was also found that compaction can be maxi- 
mized through an annealing procedure. This process in- 
volves a slow "cooling" of the system starting from a high 
shaking intensity T. These slow relaxation and annealing 
properties of this system are reminiscent of conventional 
glasses. Another qualitative similarity to glasses is ob- 
servable in the density fluctuation spectrum of the gran- 
ular system near equilibrium. The spectrum was found to 
be strongly non-Lorentzian || , revealing the existence of 
multiple time scales in the system. The shortest and the 
longest relaxation timescales differ by as much as three 
order of magnitude, and the behavior of the spectrum 
for the intermediate frequencies is highly non-trivial; in 
certain regimes it can be fitted with a power law. 

These previous experimental observations are sugges- 
tive of glassy behavior and this connection has been ex- 
plored in recent models of compaction using ideas from 
magnetic systems H. However, a more direct test of the 



glassy nature of granular compaction comes from mea- 
surements of the response of the system to sudden per- 
turbations in the effective temperature, given by T. This 
idea originates from classical experiments for the study of 
aging in glasses Q and has recently been explored using 
computer simulations ]?J . In this letter we present direct 
experimental observations of memory effects in a vibrated 
granular system obtained by measuring the short-time 
response to an instantaneous change in tapping acceler- 
ation T and propose a simple theoretical framework. 

We used the experimental set-up described in 
refs. J2|j3|]: 1 mm-diameter glass beads were vertically 
shaken in a tall, evacuated, 19 mm-diameter glass tube, 
and the packing density of the beads was measured using 
capacitors mounted at four heights along the column. 

The simplest form of this experiment consists of a sin- 
gle instantaneous change of vibration intensity from Ti to 
T 2 after t taps. For T 2 < Ti (Fig. |a) we found that on 
short time scales the compaction rate increases. This is 
in sharp contrast to what one may expect from the long- 
time behavior found in previous experiments where the 
relaxation is slower for smaller vibration accelerations. 
For r 2 > Ti (Fig. [l)b) we found that the system dilates 
immediately following tg. 

These results too, are opposite from the long-time 
behavior seen in previous experiments where the com- 
paction rate increased: Not only does the compaction 
rate decrease, it becomes negative (i.e. the system di- 
lates). Note that after several taps the "anomalous" di- 
lation ceases and there is a crossover to the "normal" 
behavior, with the relaxation rate becoming the same as 
in constant-T mode. Thus, most of the shaking history 
is forgotten after a short time. 

These data constitute a short-term memory effect: the 
future evolution of p after time t$ depends not only on 
p(to), but also information about the previous tapping 
history, contained in other "hidden" variables. In order 
to demonstrate this in a more explicit manner, we mod- 
ified the above experiment. In this second set of three 
experiments the systems was driven to the same density 
po with three different accelerations To, T\, and r 2 . After 
po was achieved at time t^, the system was tapped with 
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the same intensity Tq for all three experiments. As seen 
in Figure [|, the evolution for t > t strongly depends 
on the pre-history. The need for extra state variables in 
the problem is consistent with strongly non-Lorentzian 
behavior of the fluctuation spectrum, observed in earlier 
experiments 0. Indeed, if the evolution could be pre- 
scribed by a single master equation for local density, it 
would result in a single-relaxation-time exponential de- 
cay of the density fluctuations near equilibrium. Instead, 
a wide distribution of characteristic times is suggested by 
the spectrum. 
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FIG. 1. Evolution of the packing fraction, p, at four heights 
in the column, as a function of tap number, t. Two different 
single-switch experiments: (a) T was lowered from 5.6 to 1.8 
at to = 25 ; and (b) T was increased from 3.5 to 6.3 at to = 30. 
Curves are shifted vertically for clarity. Each curve is an 
average over 4 runs, and the measurement uncertainty in p is 
4 x 10" 4 . 



dominated by a number of local flip-flop modes with rel- 
atively high transition rates in both directions. 
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FIG. 2. The time evolution of packing fraction p for a sys- 
tem which was compacted to po = 0.613 at time to using 
three different accelerations: Fi = 1.8 (circles), To = 4.2 (tri- 
angles), and T2 = 6.3 (diamonds). After the density po was 
achieved, the system was vibrated at acceleration Fo. The 
evolution for t > to depended strongly on the pre-history. 
Each curve is an average over four experimental runs. 

This model replaces the complicated configuration 
space with a set of independent two-state systems, each 
of which is characterized by two transition rates, n e ^ g > 
I K g ^ e gives the ratio of the equilibrium prob- 
abilities of populating each state: "ground" and "ex- 
cited" . As we have argued, the higher probability ground 
state is typically the one with higher density, i.e. the vol- 
ume change v between the ground and the excited states 
is normally positive (see Fig. [| for a schematic descrip- 
tion of the model) . Our two-state approximation is close 
in its spirit to recent Grinev-Edwards and de Gennes 
models §,§. 

We now introduce the concept of a base-line density, pb, 
which corresponds to all the elementary modes at their 
ground states. Obviously, the experimentally-observed 
density is different from pi, due to a non-zero fraction of 
excited states: 



P = Pb(t) 1 



v ^ 
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To give a theoretical interpretation of the above re- 
sults, we view the problem as an evolution in the space 
of discrete "microscopic" states corresponding to differ- 
ent realizations of the packing topology (i.e. of the con- 
tact network). For each tap there is a possibility for a 
transition from one microscopic state to another. Since 
the dynamics is dissipative and the system is under ex- 
ternal gravity, a transition to a denser configuration is 
typically more probable that the reverse one. We now 
assume that the short-term dynamics of the system are 




FIG. 3. Schematic description of the two-state model: at 
every tap a transition can occur from the ground to the ex- 
cited state or vice versa. 

The summation here is performed over all the domi- 
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nant two-state modes, V is the total volume and 
the volume difference between the excited and the ground 
n-th state. Assuming that the vibration intensity T is a 
qualitative analog of temperature, we expect the popu- 
lation of the excited states, P(T) = (1 + >4"ie/ 'n^g)^ 1 , 
to grow with T, starting from zero at T = 0. Hence, for 
a given p^, the total density p will be lower at higher ac- 
celeration. This explains the observed effect of an abrupt 
change of T. After a switch from T\ to T2 at time to = 0, 
the flip-flop mode contribution to the total density would 
relax to its new equilibrium value in the following way: 



A 



rir 2 



(t) = pb / vFr lt r 2 ( v > K ) (1 — exp(— Kt)) dvdn (2) 



Here k is the relaxation rate of an individual mode, and 
the distribution function F^ 1i t 2 (v, K ) is introduced as fol- 
lows: 

Fr u r 2 (v,K) = - £ ( p(n) (r 2 ) - P (n) (ri)) 5(v - 

n 

*(«-«w e (r a )-«W fl (r 2 )). (3) 

Since the observed density changes in compaction ex- 
periment are normally less than 1% of the total density 
one can estimate the typical separation between 
neighboring flip-flop systems as 5 particle sizes, which 
yields a good justification for our no-coupling approxi- 
mation. According to Eq. (j^), if Fr 1 .r 2 does not vanish 
in the limit n — > 0, the late stage of the relaxation of 
Ar 1 r 2 {t) is given by the power law: 



r. First, we reach, using annealing dynamics, a "quasi- 
steady" state at amplitude r , so that one can consider 
Pb constant later on. Let us switch the shaking accelera- 
tion from Tq to Ti for a finite number of taps St, and then 
switch it back to Tq. During the intermediate Ti-stage, 
the system does not have enough time to completely re- 
lax to its new equilibrium. In our two-state model, the 
modes whose relaxation rate (at Ti) is below St^ 1 re- 
main unrelaxed. Assuming that the slow modes at Ti 
are mostly the same as at To, we can calculate the back- 
ward density relaxation similarly to Eq. (||), with F(v, k) 
effectively depleted below a minimal rate, kq. This cut- 
off frequency, kq, is expected to decrease monotonically 
with increasing perturbation duration St. The resulting 
density relaxation after returning to Tq is given by: 



*r 11 r a (*)sAr 1 r 9 (t)-A ri r a (oo) = 
Pb / v Fp 1 .r 2 ( v i K)e~ Kt dvdK ~ — . 



(4) 



Note that p\, is also dependent on time: although this 
cannot be described within our two-state approximation, 
the collection of elementary modes slowly evolves. Thus, 
one can observe two different processes: on short time 
scales, a fast relaxation due to the flip-flop modes is 
dominant, while over the long times, the dynamics are 
determined by the logarithmically slow evolution of the 
baseline density Pb(t). The crossover between the two 
regimes is particularly obvious in Fig. ^b, where it re- 
sults in a non-monotonic evolution. For experiments per- 
formed at sufficiently late stages of the density relaxation, 
the dynamics of the baseline density could be neglected 
compared to the contribution of the flip-flop modes (note 
that what we call a late-stage relaxation corresponds in 
fact to mesoscopic time scales which are always shorter 
than the relaxation time for pt). It has to be emphasized 
that the described experiments provide us with a tool for 
study of the response of the system, which is not limited 
to the nearly-equilibrium regime. 

One can use our simple model to predict the response 
of the system to a more complicated pattern of changes of 



3rx,r, 



,(t) 



exp(-KQt) 



(5) 



We tested these predictions by performing this three 
stage experiment, varying the duration, St, of the 
perturbation(ri) stage (Fig. ||). As predicted, the time 
needed to recover the steady-state density increases with 
the number of taps St spent in the "hot" regime T\ > Tq. 
In the coordinates chosen, the relaxation curves should 
follow the St — 00 dynamics until the saturation at 
the cut-off time, KQ 1 (St). We approximate the distri- 
bution function F by a constant above this low fre- 
quency cut-off at KQ 1 (St), up to a high-frequency cut-off, 
Kh — ltap -1 . This eliminates the unphysical low-t di- 
vergence in Eq. Figure [| shows fits of the data to 
Eq. (0), where Ko(St) is determined from the fit. The 
best-fit is achieved at Kh = 0.4, and the variation of this 
parameter would result in a simple rescaling of the time 
axis. 

Figure || demonstrates good agreement between model 
and experiment, with some systematic error at the ear- 
liest relaxation stage (which is an expected result of our 
oversimplified description of the short time dynamics). 
For the late stage relaxation, we conclude that (i) within 
our experimental precision, the St = 00 relaxation is con- 
sistent with the predicted 1 jt law; (ii) finite-5t relaxation 
curves can be parameterized by a low frequency cut-off, 
kq; and (iii) kq is a decreasing function of the waiting 
time St, shown in the insert of the Figure |j We now re- 
late our picture to previous experimental and theoretical 
results. As discussed earlier, the wide range of relaxation 
times reveals itself both in our response measurements 
and in the the fluctuation spectra of the density. It is 
tempting to relate these two kinds of data through an 
analog of a Fluctuation-Dissipation Theorem (FDT). Of 
course, there is no fundamental reason for FDT to be 
applicable to the granular system |l(J. Even though the 
above two-state model could be mapped onto a thermal 
system (in which FDT is expected to work), the ther- 
modynamic variable conjugate to density in the context 
of such a mapping has no clear physical meaning. Nev- 
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ertheless, below we outline the pseudo-FDT relationship 
expected for the granular system under rather natural ap- 
proximation. Namely, we neglect the correlation between 
the volume change v and the life time k _1 of an individ- 
ual mode, i.e. assume i<b,r(w, k) — f(n)g{v). Then the 
density autocorrelation function can be written as fol- 
lows: 

(S P (0)6p(t)) r = ^J ((v 2 ) - (v) 2 ) exp(-Kt)f( K )dK = 

( v 2\ _ ( v )2 
P^fi^-W)- («) 

Thus, the density correlator is simply proportional to 
the response function corresponding to the switch be- 
tween a very low acceleration (at which virtually all the 
modes are in their ground states) and the given one, T. 
An experimental check of this relationship requires fur- 
ther high-precision studies of both the relaxation dynam- 
ics and the fluctuation spectrum. 
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FIG. 4. The density relaxation (A = p(t) - p(ti)) of the 
system, prepared by tapping for a long time at To = 1.8 and 
then tapping for a variable number, St, of taps at a "hotter" 
intensity Fi =4.2 before being returned to To at time ti . Note 
that p(ti) is a function of St. Thus, while the asymptotic value 
of A depends on St, all curves asymptote to the same density 
p(oo), which depends only on To and not on St. The solid 
lines represent the theoretical curves, with appropriate values 
of the parameter ko. The dependence of the cut-off rate Ko 
on the waiting time St is shown at the insert for St < 4 taps. 
We do not show the value for St = 8 since we found it null 
within the error bar, as for St — 4. Each experimental graph 
is an average of 12 runs. 

Our model also gives a simple interpretation to the de- 
creasing dependence of the steady-state density on T: it 
can be attributed to the growth of the population of the 
excited states, P (T). Indeed, the corresponding correc- 
tion to the total density is about 1%, i.e. of the same 
order as the variation of the equilibrium packing fraction 
with r §. 

The slow dynamics associated with the evolution of 
the base-line density can also be addressed within our 
approach. For doing so we need to account for the 



coupling of individual modes. Namely, it would be a 
reasonable assumption that a relaxation of one mode 
to its ground state may frustrate such a transition for 
some of its neighbors (e.g. in 3D the most compact lo- 
cal cluster can be created only at the expense of less 
dense neighboring regions). Thus, we arrive at an ef- 
fective anti-ferromagnetic (AF) coupling (of an infinite 
strength) between the two-state modes. This extension 
of our model makes it remarkably similar to the so-called 
reversible Parking Lot Model (PLM) 0JT§], which has 
been successful in describing many aspects of granular 
compaction experiments PJjl. Recent simulations based 
on the "tetris model" Q for compaction also find slow 
glassy responses to changes in T, but do not capture the 
short-term memory effect described here. p~4| . 

In the PLM, D-dimensional space (the parking lot) gets 
packed with finite-size objects (cars), which may arrive 
and depart with fixed rates and which are not allowed to 
overlap. Now, the transcendental relationship between 
PLM and the granular compaction experiment is easier 
to explain: both PLM and our coupled flip-flop model 
belong to the same generic class of AF-type systems (in 
the case of PLM, a local two-state mode is represented 
by a particle whose center of mass may or may not be 
placed at point x; the mode coupling is due to the hard- 
core interactions). The PLM is known to capture the 
slow dynamics of granular compaction and some features 
of its fluctuation spectrum PfT^]. In fact, we performed 
numerical simulations of the PLM that display [|l3| the 
memory effects discussed in this work. 

In conclusion, we used a sequence of abrupt switches 
of the shaking intensity T to study the response of a vi- 
brated granular system. This technique can be used in 
the vicinity of the steady state density, as well as far 
from equilibrium. The major result is the direct demon- 
stration of a memory effect in the system: the evolution 
is not predetermined by the local density alone, and its 
description requires introduction of additional "hidden" 
variables. Our phcnomcnological model for this behavior 
is built on minimal assumptions about the dynamics of 
the system and produces results which are generic and 
are expected to be valid for a wide class of more realistic 
microscopic models. 
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